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We describe in detail a physical situation in which instantons are necessarily complex, not just 
Wick rotations of classical solutions to Euclidean spacetime. These complex instantons arise in the 
semiclassical evaluation of vacuum pair production rates, based on Feynman's worldline path integral 
formulation. Even though the path integral is a sum over all real closed trajectories in spacetime, 
the semiclassical description of non-perturbative pair production is dominated by closed classical 
trajectories that are generically complex. These closed trajectories contain segments associated 
with nonperturbative instanton suppression factors as well as segments producing phase factors that 
incorporate quantum interference effects. For a class of time-dependent electric fields we implement 
this procedure and demonstrate excellent quantitative agreement with alternative methods. 

PACS numbers: ll.15.Kc, 12.20.Ds, ll.15.Tk, 03.65.Sq, 

I. INTRODUCTION 

The Heisenberg-Schwinger effect is the non-perturbative production of electron-positron pairs from the quantum 
electrodynamical (QED) vacuum under the influence of an external electric field [iHS]. The general quantum field 
theoretic formalism for computing the associated probability was developed by Schwinger in terms of the QED effective 
action [3]. However, there are still serious obstacles to the implementation of a reliable quantitative computation when 
the external electric field is taken to be that for realistic high-intensity laser pulses. Interest in this problem has been 
revived recently, spurred by new experimental developments in ultra- high intensity lasers [6] . Models of laser pulses 
with one-dimensional inhomogeneities, such as time-dependent linearly polarized electric fields are well understood 
(although the question of pulse sequence optimization still stands), but the situation is much less clear for fields with 
multi-dimensional inhomogeneities, such as occur naturally in more realistic physical configurations of colliding high- 
intensity, spatially focussed, laser pulses [7j. This is a pressing matter, because recent theoretical progress suggests 
[5] that the critical peak field intensity required to observe this effect may in fact be several orders of magnitude 
lower than the estimate based on assuming a constant electric field [TJ [2], raising hopes that the effect may be 
observed experimentally in the not too distant future. In turn, this also raises important unresolved questions about 
back-reaction and cascading effects [S]. 

In the quantum field theoretic approach [2, the technical problem is to compute the non-perturbative imaginary 
part of the "effective action", r[A] — —ih In det [mc — ip], where the Dirac operator, p= ^^{d^ — i-^A^j), defines the 
coupling between electrons and the applied (classical) electromagnetic field that represents the field produced by 
the laser pulse. The conventional approach to this problem in the case of a one-dimensional inhomogeneity reduces 
it to a Id scattering problem |10nl2j . invoking Feynman's picture of anti-particles as particles traveling backward in 
time [IST. There are then many possible approaches to compute pair production rates and the momentum spectra of 
the produced particles [T0Hl2l[T4Ul7j . However, these one dimensional methods do not generalize in a simple, efficient 
way to the multidimensional situation. There have been recent developments for multidimensional fields concerning 
finite-plane- wave fields [TB], and the numerical implementation [T5] of the Dirac-Heisenberg-Wigner formalism |20j . 
On the other hand, a natural semiclassical formulation of the general problem is in terms of worldline instantons, 
a semiclassical approximation to Feynman's worldline path integral expression for the QED effective action. This 
method has been quantitatively confirmed for certain one-dimensional field configurations, and the general formalism 
has been outlined for multi-dimensional field configurations pT| - B5] . A technical obstacle to the implementation of 
the worldline instanton method in higher dimensions has been the physical interpretation of the complex classical 
trajectories that naturally arise. The purpose of this paper is to clarify the physical meaning of such complex classical 
trajectories, using a one-dimensional example for which we can confirm our results by comparison with other methods. 

Usually instantons appear as solutions to the Euclidean classical equations of motion, in which ~^ — ix^. In 
fact, this definition is too restrictive for the worldline picture, and a more natural definition is to seek solutions with 
imaginary proper-time: r — > s = ir, as proposed by Rubakov et al |24) . In simple text-book cases this transformation 
to imaginary proper-time goes hand-in-hand with the Wick rotation to imaginary (Euclidean) time, but there are 
examples in which the spacetime instanton trajectories x^(t) should be viewed as lying in complex Minkowski space 
[211 [25]. In the case of QED, as studied here, the situation is even more interesting because the gauge coupling 
produces a Lorentz-force term in the relativistic classical equations of motion, a;^ = F^i,{x)x'^ , which acquires a 
factor of "i" after rotating to imaginary proper time, so that the instanton equations are manifestly complex from 
the very beginning. [This is analogous to the effect of a magnetic field on a tunneling problem in non-relativistic 
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quantum mechanics |26| : it breaks time-reversal symmetry and makes the tunnehng instanton equations complex.] 
In this paper we show that for the problem of QED vacuum pair production, complex instantons are needed to 
capture the physics of quantum interference between distinct instanton trajectories. This phenomenon of quantum 
interference arises for laser pulses with temporally localized electric field pulse shapes having sub-cycle structure, 
such as "carrier-envelope-phase" or "chirp" features [?7H2^ . In addition, appropriately chosen temporal sequences of 
pulses can produce significant coherent enhancement in certain momentum modes, an explicit time-domain realization 
of multiple-slit interference . We treat both scalar and spinor QED to show explicitly how the interference terms 
are affected by the quantum statistics of the particles. We also note as motivation for studying complex instantons in 
QED the fact that complex trajectories are well-known in multi-dimensional tunneling phenomena in non-relativistic 
quantum mechanics |31H35j . Furthermore, the physical meaning of complex classical trajectories has recently been 
further elucidated by the study of PT-symmetry in quantum mechanics 36J . 

In Section II we recall the worldline instanton formalism for the QED effective action, and explain why complex 
instanton solutions appear. In Section III we present the worldline instanton solution for the more general problem 
of finding the momentum spectrum for the produced electron-positron pair, and state the appropriate boundary 
conditions for finding the semiclassically important solutions of the complex classical equations of motion. Quantitative 
results are presented in Section IV, demonstrating excellent agreement with alternative methods of solution, and 
Section V contains our conclusions. An Appendix discusses an important and interesting numerical instability that 
occurs for certain ranges of values of the longitudinal momentum, and we present a simple resolution of this instability 
by taking advantage of the reparametrization invariance of the worldline path integral. 



II. WORDLINE INSTANTON FORMALISM 

A. Worldline form of the QED effective action 

Following Schwinger [3], we compute the non-perturbative pair production probability P from the imaginary part 
of the effective action rcff[^] for the QED vacuum in a prescribed classical background field A^t(a;): 

P=l- e-2im[roff]/n ^ 2 j^^^^ 

For physically relevant configurations, ImFoff [A] /h is extremely small, which justifies the approximation in the last 
step. The effective action Foff [A] is defined, for spinor and scalar QED respectively, as (henceforth we set h = 1) [5]: 

F^P'"°'- [A] ^ ~i In det [mc ~ip]^~^tr\n [m^c^ + p^] (2) 
pscaiar = I In dct [m^c^ + = I tr In [m^c^ -|- (3) 

The covariant derivative operator 13^ has been defined above, in the Introduction, and we adopt the space-time 
metric convention = diag(l, —1, —1, —1). Both Schwinger [3] and Feynman [37] interpreted these effective actions 
in terms of quantum mechanical propagation in four-dimensional spacetime: 



r:r-[A] = e-^^tre-«-^ , H., ^ (4) 







The factor of 1/2 in "H is a convention |37j . introduced by simple analogy with the form of the Hamiltonian in non- 
relativistic quantum mechanics. The integration variable T can be thought of as the total propagation "time", which 
leads naturally [37] to a path integral expression for the effective action. For scalar QED: 

pscalarj^j = -z H ^g-^^^T f ^4^ f Vx C^'^^^^ (6) 

Jo ^ J Jx{0)=x(T) 

where S is the classical action for a relativistic scalar charged particle, coupled to the gauge field Af^{x), propagating 
around the closed trajectory x^{u) with a propagation period T: 
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The paths x^{u) are closed paths in four-dimensional spacetime, parametrized by a parameter u, which we relate 
to proper-time in the following subsection. For spinor QED there is an additional spin interaction, because = 
D^— ^(J^'^F^jj, and the effective action for spinor QED can be written with an additional Grassmann path integration 

2 Jo ^ J Jx(0)=x{T) J 

where 5 is as in ([t]). 



B. Semiclassical approximation: worldline instantons 

The path integrals in ([6| and Q are of course only analytically calculable in special idealized cases, so we must 
resort to approximation methods. For non-pertubative questions such as the the pair production probability one can 
use the numerical worldline approach [41 , or a semiclassical evaluation based on worldline instantons [HI [22]. In this 
paper we follow the worldline instanton approach, in which we search first for a saddle point solution to the bosonic 
path integral by solving the classical equations of motion for relativistic motion of a charged spinless particle: 

^ = (9) 
dw^ c du 

A closed trajectory solution to these classical equations is called a "worldline instanton" . In certain physical situations, 
such as for non-perturbative pair production from vacuum, these classical solutions give a dominant contribution to 
the path integral in ^ and_p]). 

The classical equations l\9« have an obvious first integral, the "energy" H — p^,^ — L — 5-3^^^, which is a 
constant of motion. This constant is fixed by making also a saddle-point approximation to the T integral, which gives 
a critical condition: 

m^c^ dS „ 

Since the variation of the action with respect to the period T, namely is equal to minus the conserved energy, 
this implies the normalization 

dXfj^ dx^ 22 /1 1 \ 

Using the relation, ^ = 1/ \/l — t/^/c^, between time t = x^/c and proper-time r, this identifies the propagation 
parameter u as a multiple of the proper-time: 

u=- (12) 

m 

We can therefore identify raT with the total proper-time of the quantum mechanical evolution in ([5]). [We shall 
see later, in the Appendix, that a different scaling of the proper-time evolution parameter leads to some numerical 
advantages in certain situations.] 



The critical saddle point period Tc is determined by the condition (10), and when we evaluate the full exponent of 
the T integral in ([g]) we obtain Hamilton's characteristic function W 



with a prescribed "energy", equal to ^m?c^. This is just the familiar Legendre transform of classical mechanics, 
relating the action S[x]T] and the characteristic function = S[x]T] + £T, which implies: = and 



refer to this as the classical action for the motion of the relativistic charged particle to traverse its closed trajectory 

dT 

= T. Here, in this relativistic problem the role of the "energy" £ is played by |to^c^, in complete agreement with 
Feynman's interpretation of the Klein-Gordon and Dirac equation in terms of proper-time |37j (based on results of 
Fock and Nambu [42l[43]). Thus, on the classical solution, the characteristic function becomes 

, ,, , , 1 no, /""^ / 1 dx,, dx^ e dx,, , ,, , , ^ o o\ dx^ 

^ ^ ^' 2 ^ n \2 du du c du ^ ' 2 J L du ^ ^ 



This is the classical function that appears in the exponent after making a semiclassical approximation. Note that 
there may be several classical saddle point trajectories x^{u) relevant to the semiclassical approximation. 
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C. The need for complex worldline instantons 

The worldline instanton approach was originally suggested for QED with a constant background electric field |21| . 
then extended to QED in background fields with one-dimensional inhomogeneities in j22j , and a general formalism was 
also proposed for more general fields, based on an analogy with Gutzwiller's trace formula The most difficult part 
of the computation is to find the semiclassically important classical trajectories. The first observation is that we seek 
closed trajectory solutions; this is because the effective action involves a trace, so that the quantum mechanical path 
integral in ([6j [S]) is expressed as a sum over closed paths in four-dimensional spacetime. But we still need to specify 
some initial conditions in order to search for appropriate saddle point solutions. Previously it had been suggested 
to look for closed trajectory solutions to the Euclidean classical equations of motion [521 [53], obtained from ([9]) by a 
Wick rotation, ^ x'^ = ix^. We show here that this prescription needs to be refined and extended, in order to 
describe quantum interference effects. 

We begin with the formalism of Rubakov et al [Ml[3Sll31j, that an instanton solution is associated with a deformation 
of the contour of the T integral onto the imaginary axis, so that we look for solutions with imaginary "proper-time" 
parameter: u s — iu. In the familiar case of scalar [non-derivative] couplings, the classical equations of motion 
acquire a sign change: 

d^a;^ dV(x) d'^x>' dVix) 

^ = (14) 



dv? dx^ ds-^ dx^ 

Note that the equations of motion remain real under this operation. On the other hand, for the gauge coupling of the 
QED case, imaginary proper-time introduces a factor of " i" into the classical equations of motion. 

^ = ^ ^ i^- F^^^x)'^ (15) 

atf^ c du ds-^ c ds 

Therefore, the instanton equations are generically complex, and so the solutions will be generically complex, as will 
be the classical action evaluated on such a solution. This raises the question: what is the physical significance of such 
complex classical solutions? We answer this question in the remainder of this paper. 



FIG. 1: Sketch of basic shapes of gauge fields [solid, blue curves] A3 (a;") and the corresponding electric field [dotted, red, curves] 
Es{x°). In the first case, the gauge field is a monotonic odd function of time and the electric field is a single localized pulse. 
In this case there are no quantum interference effects. For the next two cases, the gauge field is an even function of time, the 
electric field is an odd function of time, and there are significant quantum interference effects, as discussed in [28j using WKB. 

An important comment concerning previous work on worldline instantons |21l I22j , as well as on early work on the 
WKB scattering approach [T0Hl2] . is that they concentrated on two particular sub-classes of problems. The first class 
consists of time-dependent electric fields linearly polarized along a particular spatial axis [say, the a;^ axis] , so that the 
gauge field can be written with only one non-zero component ^3(2;''). Classic cases are: ^3(0;'^) = Ex'^, for a uniform 
field; ^3(0:°) = E /uj sm{ujx'^) for a monochromatic sinusoidal field; A3(j:°) — E /uj ta,nh{ijjx'^) for a temporally localized 
single-pulse field. All these examples have the important property that v43(x°) is an odd function of time. The second 
class of fields involves static but spatially inhomogeneous fields represented by a scalar potential Aq{x). Now, observe 



that in each of these cases, the complex classical equations of motion in ( 15 ) reduce again to real equations if we 
combine the analytic continuation, u s — iu, with the Wick rotation, a;° — > = ix^, and A'^ — J- A* = iA'^. In the 
former case, this is only true if A^i^x*^) is an odd function. Otherwise, the equations remain complex. This explains 
why the previous analyses were able to produce consistent and correct results using as worldline instanton equations 
the classical Euclidean equations of motion with imaginary proper-time. 

As a simple illustrative example, consider the case of a constant electric field, of strength E. Then the classical 
trajectories are hyperbolic paths: x^{u) = sinh{eEu/c), and x^{u) = cosh{eEu/c). These are not periodic 
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closed paths, but we can make them closed and periodic if we take u ^ s = iu. We then obtain complex solutions, 
= —i sm{eEs/c), and = cos{eEs/c), which become real again when expressed in terms of the Euclidean 
time = ix^. But the reality of the worldline instanton solution in Euclidean spacetime is an "accident", a direct 
result of ^3(2;°) being an odd function of x'^ . For this solution, the period is = and we can evaluate the 



characteristic function ( 13 1 as: 



(16) 



Then the semiclassical approximation to the scalar QED effective action leads to 

Im r^'jf i'^'' K V e-'^'-'"'™ KVe- (17) 

which is the familiar result of Sauter ^j, Euler and Heisenberg [2] and Schwinger [3], and is a simple prefactor. 
Other examples in which the explicit worldline instanton trajectories and actions can be evaluated in closed form are 
given in ^TD , with results in agreement with the corresponding WKB treatment [rC.TT] (including also the prefactors). 

In more realistic time-dependent electric fields, such as those having an envelope structure as well as an oscillatory 
structure, there are quantum interference effects, which can produce both enhancement and suppression |27j . These 
cases are associated with vector potentials that cannot be written as an odd function of time, as sketched in Figure 
[ij The semiclassical analysis of such systems, incorporating quantum interference, has been given in |28j using the 
WKB approach. Here we explain how to solve these quantum interference problems using the more general formalism 
of worldline instantons. Since the worldline instanton equations remain manifestly complex, we confront directly the 
problem of the physical meaning of complex instantons. Our motivation for studying this class of problems first is 
that we have results with which to compare, so that we can quantitatively verify the validity of our approach. 

There may be several classical saddle point trajectories, x^--^{u) , labeled by an index (j), and the imaginary part 
of the effective action is then approximated by 

Im r^^^'^' « r'-^^ '^in.^t.nton (ig) 

j 

where W^^^^^^^^^^^ is the characteristic function (13), evaluated on the j"^ saddle point solution, and V^^'' is a (readily 
calculable) prefactor. In the case of spinor QEDTThe saddle point trajectory is the same, but there is an additional spin 
factor coming from the evaluation of the (Gaussian) spinor path integral in ([s]), evaluated on the critical trajectory 

x'-^iu) m-- 

Imr~«-^pWdetV^ e-^"^— (19) 

The determinant spin factor can be computed straightforwardly using the Gelfand-Yaglom method, since it only 
involves ordinary differential operators, as can the prefactors HO] , 



III. MOMENTUM SPECTRA FOR PARTICLES PRODUCED IN TIME-DEPENDENT ELECTRIC 

FIELDS 



A. General formalism 



In this section we extend the worldline instanton method of f5TJ [51] to compute not only the total probability 
of pair production, but also the momentum spectrum of the produced particles. In the situation where the vector 
potential A^ixP) is a function only of x", the spatial momenta of the electron-positron pair, p = (PijPs), are good 
quantum numbers and can be used to characterize the final states. Furthermore, in a strong field the pair production 
is predominantly along the direction of the electric field, so we can neglect p_L and concentrate on the dependence 
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of the number of produced pairs on the longitudinal momentum p^. Studies of temporally structured electric field 
pulses have revealed an intricate dependence on pa, due to quantum interference effects [27J[28]. In this Section we 
show that the analysis of this momentum dependence using worldline instantons requires complex worldline instanton 
trajectories, not simply Euclidean classical trajectories. 

To address the momentum dependence we convert the worldline path integral expressions ([6j [s]) for the QED 
effective action into phase space path integral expressions. Since the only difference between spinor and scalar QED 
in the semiclassical approximation is the spin factor determinant in (19), we first concentrate on the scalar QED case 
in order to find the instanton trajectories. 



pscalar f 



dT 







x{0)=x{T) 



off '[A 

where the classical Hamiltonian density is 

H{x,p) = p^'^ ~ L{x,i) ^^(^Pf, + ^A^{x)^ 



(20) 



(21) 



Since the external field yl^ is independent of the spatial coordinates x, the spatial path integral can be done, producing 
delta functions in ^ , thus imposing the conservation of spatial momentum. This means that the functional integrals 
for the spatial momentum reduce to ordinary integrals: J D^p — J d^p. Then we can convert the remaining phase 
space path integral over x° and p^ back to a coordinate space path integral, leading to a worldline path integral 
expression in terms of a single coordinate (u) , parametrized by the spatial momenta p: 



'[A] ^-ilV^ I dW I dp3 ] r 



T 



^ I dx'^ 2?a;°e-^^[-"l 

Jx°{0)=xO{T) 



where the classical action is: 



(22) 



(23) 



Recalling the scaling (12) between between u and proper-time t, this effective action expression (22 23 1 has the form 
of a quantum mechanical path integral in the single (time) coordinate x^(t), parametrized by the proper-time r, with 

a "potential", V(x°) = —\ (pa -|- |A3(x°))^ — that depends parametrically on the spatial momenta p^ and pi_. 
This agrees completely with the WKB picture of pair production for time dependent electric fields as a quantum 
mechanical Schrodinger over-the-barrier scattering problem in the time-domain [101 lll[ I28j . The argument of the 
p integral in (22 1 has a non-perturbative imaginary part that gives the probability of producing pairs with spatial 
momentum p. 



B. Semiclassical approximation for momentum spectrum 



We now make a semiclassical approximation to the effective action in (22 23). First, we neglect the transverse 



momenta, setting p±^ = 0. This only affects prefactor terms, and can easily be incorporated if desired, and will not 
be important in what follows. Then the classical equation of motion is 



(24) 



This equation (24) is of course just the remaining nontrivial classical equation of motion from (|9]) after integrating 
the x^ equation to give = Pz + f ^3(a^°), with p^ arising as an integration constant. In [21], in computing the total 
pair production rate rather than the momentum dependence, this integration constant was taken to vanish, with a 
Gaussian momentum integration producing certain prefactor terms that were explicitly computed and shown to agree 
with the WKB result [TT]. Here we retain the p^ dependence in the equation of motion (24) in order to find the 
longitudinal momentum spectrum of the produced particles. Thus, the solution x'^{u) will depend parametrically on 
P3, as will the classical action when evaluated on that classical solution: 



S[x\u);T] = 




rn^c^ I du 

2 I 



(25) 
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Here we have set p± = and used the proper-time relation 

\ c / 



(26) 



which also expresses the existence of a first integral for the equation of motion (24). The term —^m?c^T in (25) 



cancels against a similar term in the T integral in ( 22 ) , and so the resulting exponent is the Hamilton characteristic 

1 

2 



function, the action for a closed trajectory of fixed "energy" ^m^c^: 



Then the imaginary part of the QED effective action, in scalar and spinor QED respectively, is given by: 



(27) 



(28) 



Imr~ « - / dpsE^^'^det^/^ fs,^l^-^^^F,,i4\u;p,))] e^^/— t^^) 
where the sum is over all relevant semiclassical trajectories (to be specified explicitly in the next subsection). 



(29) 



C. Boundary conditions for the worldline instanton trajectories 



We now specify the appropriate boundary conditions for finding solutions to the classical equations of motion ( 24 1 



To find instanton solutions we expect to take u (and hence T) to be imaginary. This corresponds to a deformation 
of the contour of integration of the T integral [M]. This then leads to (some) closed path trajectories. However, for 
background gauge fields with nontrivial temporal structure, there can be more than one different instanton solution, 
and we need to be able to find all relevant instanton trajectories and tie them together. Thus, we consider a contour 
for u that can have "vertical" segments along which the real part is constant, which we refer to as instanton segments, 
and in addition we consider "horizontal" segments along which the imaginary part of u is constant, which we refer to 
as "interference" segments. This situation is illustrated in Figure [2] This leads to a difference of sign in the equation 
of motion ( 24 ) for the two types of solution. 



|Im[u] 



(1) 



(3) 



~1 



Re[u] 



(2) 



FIG. 2: Sketch of the complex contour in the complex u plane, showing two distinct "instanton segments", labeled (1) and (2), 
having constant real part of u, and an "interference segment", labeled (3), having constant imaginary part of u, and connecting 
the two distinct instanton segments. 

Next we specify the initial conditions. As mentioned at the end of Section HIA, the momentum spectrum problem is 
equivalent to the WKB scattering problem treated in [23". Motivated by this WKB analysis, we propose the boundary 
condition that the classical worldline trajectories should begin and end at WKB turning points, which are defined as 
points where 



mV+(p3 + ^A3(xO))' = 



(30) 
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These points lie in the complex x° plane, and they occur in complex conjugate pairs in the physically relevant case 
where A3 is a real function of its argument. The turning points move around in the complex pla ne as ps varies, but 
always remain in complex conjugate pairs. Given this initial condition for a;*^, then because of (26), the corresponding 

initial condition is that 



must vanish. Thus we are led to the initial condition for our semiclassical trajectories: 



that the initial "velocity" vanishes, 
solutions as follows: 



du 



= 0, and the initial (and final) point is a turning point. We therefore seek 



• Instanton segments: For an instanton segment, we take a complex proper-time evolution parameter u with 
constant real part. Re (it) = constant, corresponding to one of the "vertical" segments in the contour depicted 
in Figure [2j Therefore, we have the following equations of motion to solve: 



a;"(Uinitial) 

du 



(p3 + ^A3(a;"))aoA3(a:°) 



a turning point solution of ( 30 ) 



(31) 
(32) 

(33) 



• Interference segments: For an interference segment, we take a complex proper-time parameter u with con- 
stant imaginary part, Im(u) — constant, corresponding to one of the "horizontal" segments in the contour 
depicted in Figure [2j Therefore, we have the following equations of motion to solve: 



du'^ 

a;"(uinitial) 

dx° 
du 



= +^fp, + ^A,ix'^))doA,ix°) 
c \ c / 



a turning point solution of ( 30 ) 



= 



(34) 
(35) 

(36) 



Notice the different sign in the equations of motion for the two types of segment. 



ImCx-^O) 



ReCx-^O) 




FIG. 3: Sketch of basic shapes of worldline instanton trajectory segments. The blue segments are instanton segments and the 
red ones are interference segments. In the first plot, we show the turning point pair for two difi^erent values of pa. When ps = 
the pair lie on the imaginary axis, and the instanton trajectory lies entirely on the imaginary axis. This is the usual 
Wick rotation to imaginary time. The other curve shows the situation for pa 7^ 0, in which case the turning points lie ofi^ the 
imaginary axis [note they still form a complex conjugate pair], and the instanton trajectory connecting them is curved. So this 
solution is not captured by a Wick rotation. In the second plot, we show the situation for a field with four turning points, 
in two complex conjugate pairs. Here, even for pa = 0, the turning points lie off the imaginary a^xis, and there are two types 
of trajectories. The blue instanton trajectory connects a complex conjugate pair of turning points, while the red interference 
trajectory connects two different turning points, not a complex conjugate pair. 



The important observation is that these trajectories go from one turning point to another. For an instanton segment 
the trajectory goes from one turning point to its complex conjugate, while for an interference segment the trajectory 
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goes from one turning point to another distinct turning point with different real part [if the field A^{x^) is such that 
there is another such distinct turning point] . This is illustrated in Figure |3] In Figure [3^ we see two different types 
of instanton segments for the case ^13(0;'') = E /oj ta,nh{ujx'^). When = 0, the turning points lie on the imaginary 
x° axis, and the instanton trajectory goes along the imaginary axis, connecting the turning point to its complex 
conjugate. This explains why the Wick rotation to imaginary time, — ix^, is sufficient for this case ,22) . However, 
when p3 0, two things change: first, the turning points move off the imaginary axis into the complex plane, and 
second, the trajectory is no longer linear. Thus, in this case with p3 7^ we must consider truly complex instanton 
trajectories for x'^{s), even though for this vector potential there are no interference trajectories (reflecting the simple 
single-bump structure of the corresponding electric field). 

The second example, in Figure pp, is a case in which interference effects do arise, taking the example ^3(2;°) = 
E/uj/ (1 + (wa:°)^), that was studied in [55] using WKB methods. Here, there are four turning points, in two complex 
conjugate pairs, and they never all lie on the imaginary x'^ axis. Thus, the naive Wick rotation to imaginary time 
is not sufficient for this problem. Nevertheless, using the prescription described above we find both instanton and 
interference trajectories, as illustrated in Figure [3|d. There are instanton segments connecting complex conjugate 
turning points, and interference trajectories connecting the two distinct turning points having different real parts. As 
P3 changes, all four turning points move, and the shape of the trajectories change, but the pattern remains the same. 





FIG. 4: Sketch of the four different complex instantons for the field in (41 1, which has four complex turning points in the 
complex x° plane, shown as solid circles. In the first two figures, the instanton goes from a turning point to its complex 
conjugate and back again, while in the last two figures the closed trajectory is composed of both instanton and interference 
trajectories. The difference between the last two is the sense of propagation around the loop. 



The total contribution to the imaginary part of the effective action is a sum of various types of closed trajectories, 
composed of both instanton and interference segments. It is best illustrated by the example in Figure |4j again for 
the field ^3(2;°) = E/uj/^l + (wa;^)^), which has two pairs of complex conjugate turning points. We can form closed 
instanton trajectories by going from one turning point to its complex conjugate and then back again. This is the 
usual textbook [33] instanton contribution e~l'^'='l to ImF, corresponding to the first two plots in Figurel4l There 
are two such contributions because there are two such instantons, one for each complex conjugate pair. [It is an 
(intentional) accident of this example that the two instantons give exactly the same contribution - for other fields 

with two separate instantons, one may have a dominant e"'*^":!'' factor.]. But there are other closed trajectories 
composed of both instanton and interference terms, forming a closed loop through all four turning points, as shown 
in the last two plots in Figure |4] . The sum over closed trajectories in (28 29) corresponds in this case to the sum 
over the four distinct closed trajectories shown in Figure [4] The generalization to fields with more pairs of turning 
points is clear, and follows the same pattern as the WKB approach in [28] . Note that closed trajectories consisting of 
an interference trajectory and its reverse do not contribute to the imaginary part of the effective action. 

For spinor QED, the only change is that we must include the determinant factors appearing in (29). For this class 



of fields, these determinant factors reduce to factors of ±1. This is because we can write the spin factor as 



tr 



4 cos 




doA3{xQ)du 



(37) 



The integral around the closed loop can be separated into instanton and interference segments, and one finds that 

(38) 



doA3{xo)du = - 



For an instanton segment this follows from a substitution y = i (ps + ^A3(x°)), while for an interference segment we 
substitute y = (p^ + ^^3(2;°)). This gives a net result of cos(m7r/2) for a closed trajectory with m segments. Thus, 
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for example in the first two trajectories of Figure [4] there are two segments, so there is a spin factor minus sign, which 
combines with an overall factor of (—1) to give a positive imaginary part of the spinor effective action. On the other 
hand, for the latter two trajectories shown in Figure |4] there are four segments, so there is a positive spin factor, and 
so these trajectories contribute with the opposite sign. For closed trajectories containing 2m segments we get a net 
spin factor of cos(m7r). This, along with the global sign difference between scalar and spinor QED is a reflection of 
the role of quantum statistics in determining the sign of quantum interference terms, and the pattern of signs agrees 
with the pattern found in the WKB treatment of quantum interference in vacuum pair production |28) . 



IV. QUANTITATIVE RESULTS 



A. One dominant pair of turning points 



In this section we compute the momentum spectrum of pairs created by a potential with just one dominant pair of 
turning points. We choose the common example of an exactly soluble gauge field: 



A{t) 



Eo 



tanh (ujt) 



E{t) = Eo sech^ (wi) 



(39) 



While this gauge field yields an infinite number of turning points in the complex plane, the semiclassical amplitude 
is dominated by one single pair of turning points lying closest to the real axis [28J . There is therefore no quantum 
interference, and neglecting prefactors, the expectation value of the number of pairs produced with longitudinal 
momentum follows from (28 29 ) as 



(40) 



where Winstanton is evaluated using (27), and iVFinstanton(P3) is real and positive. Whenpa = there are simple explicit 
formulae for the the closed trajectories [22J, but when ps the complex classical closed trajectories are given in 
terms of inverse elliptic functions, which are cumbersome. Instead, we use a direct numerical integration, for each 
P3, of the classical equations of motion ( [sT] [32| [33| . A typical trajectory is shown in Figure [s] We then compute the 
classical action Winstanton in ( 27 ) , evaluated on this classical solution. The resulting approximate expression ( 40 ) for 
the particle number is plotted in Figure [6] as a function of the longitudinal momentum p^, showing excellent agreement 
with other methods (as discussed in [28 ). 



B. Two pairs of turning points 



A more interesting example is provided by a field that induces quantum interference amongst the produced particles. 
As in [25] , we consider the field 



A{t) 



Eo 



1 



E{t) 



2EoLut 



(41) 



which has precisely two pairs of (complex conjugate) turning points in the complex x'^ plane. For each value of 
longitudinal momentum p3, we integrate the classical equations of motion with their initial conditions, as given in 
both (31 32 33) and in ([34j |35j |36]). This procedure produces four different types of closed trajectories, of the 
form sketched in Figure [4] A numerical example of an interference trajectory is shown in Figure [T] For each such 
trajectory, we then compute the classical action VFinstanton in (27), evaluated on this classical solution. Asps varies, the 



turning points move, so each trajectory also changes, as does the classical action. For the instanton-type trajectories, 
^Winstanton is real and positive, while for the interference trajectories iWinstanton has both real and imaginary parts, 
coming from the instanton and interference segments respectively. From ( 28 ) and ( 29 ) we obtain the the semiclassical 
approximations: 



(3) 



-i W 



„(P3) 



(42) 

AA"P'""(P3) « e-'^.-t-to„(P3) + e-'^.-t-"t°"(f^) - e-'^.-t-to„(P3) _ g"* >^.nia„t„„(P3) (43) 

where the superscripts label the trajectory type as shown in Figure |4] Note that for the spinor case the determinant 
factors in (29 1 are +1 for the purely instanton closed trajectories, but are equal to -1 for the interference trajectories. 
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Im(x3) 





FIG. 5: The first plot shows the complex instanton trajectory in the complex x'^ plane for the single-pulse electric field in 
(39 1. The trajectory goes from one turning point to its complex conjugate and then back again. The second plot shows a 3D 
representation of the same trajectory. The field parameters are: Eq — 0.1, uj = 0.1, and pa — 0.5, all in units set by the electron 
mass scale m. 



N 




"To Pi 



FIG. 6: The expected particle number, as a function of longitudinal momentum ps, for both scalar and spinor QED. The 
dashed (red) line shows result ( |40[ ) evaluated on the complex worldline instanton trajectory, while and thick (blue) line shows 
the WKB result from [28|. The field parameters are: Eq = 0.1, and lu = 0.1, in units set by the electron mass scale m. 



as explained in the previous section and also confirmed numerically. This encodes the difference between spinor and 
scalar QED in the semiclassical worldline instanton approximation. 

In fact, for this particular example, the different instanton segments are symmetrical, so that iW^^"^ — iW'-^'^ = 
Re(jW^(3)) ^ Re(iW/(4))^ and lm{iW'-^^) = -Im(il^(4))^ Therefore, in this case we can write 



AA--'-(p3) « 4e-<^--„(P3) eos2 (^Im (^^ ^'i^^onfe) 



instanton 



(44) 
(45) 



These expressions are plotted in Figure [8j as functions of the longitudinal momentum p^, showing excellent agreement 
with the exact numerical result and with the WKB results from 1281. 
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FIG. 7: The first plot shows a closed trajectory, in the complex x plane, with both instanton and interference segments. The 
second plot shows a 3D representation of the same trajectory. The field parameters are: Eq = 0.1, uj = 0.1, and pa = 1.85, all 
in units set by the electron mass scale m. 



In the course of this analysis we found an interesting numerical instability that arises for certain values of p3, 
when the classical trajectories approach poles in the complex plane [which typically arise for localized fields]. In the 
appendix we present a simple modification to the numerical procedure that avoids this instability, taking advantage of 
the einbein formulation of the worldline effective action, due to reparameterization invariance of the associated path 
integral. 
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FIG. 8: The expected particle number, as a function of longitudinal momentum pa, for both scalar and spinor QED. The dot- 
dashed (black) line shows the result of an exact numerical computation. The solid (red) line shows the result ( 40 1 evaluated on 
the complex worldline instanton trajectory, while and dashed (blue) line shows the WKB result from [28]. The field parameters 
are: Eq = 0.1, and uj = 0.1, in units set by the electron mass scale m. The first plot is for spinor QED, and the second for 
scalar QED. 
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V. CONCLUSIONS 

To conclude, we have shown that it is necessary to consider complex classical trajectories in the semiclassical 
approximation of the worldline path integral expression for the QED effective action, in order to obtain the pair 
production probability and the associated momentum spectrum of the produced pairs. This is somewhat surprising, 
because the original path integral is of course a sum over real closed trajectories in Minkowski spacetime. To extract 
the momentum dependence and to incorporate the effects of quantum interference it is not sufficient to consider 
classical solutions in Euclidean spacetime (where time is pure imaginary); one must consider genuinely complex 
classical solutions. The semiclassically relevant closed classical trajectories, the "worldline instantons", are composed 
of segments that we have called "instanton segments", for which the proper-time parameter has constant real part, 
and "interference segments", for which the proper-time parameter has constant imaginary part. Generically, all these 
closed trajectories are complex, and so do not arise as solutions to the Euclidean classical equations of motion. With 
just a single instanton segment this is precisely the instanton prescription of Rubakov et al |24j . If there is more than 
one instanton solution, then in addition we need to include the interference segments that tie the instanton segments 
together. These interference trajectories, as shown in the last two plots in Figure |4] produce classical actions iW that 
have both real and imaginary parts, with the imaginary parts characterizing the quantum interference effects. For 
spinor QED there are additional determinant factors that give additional minus signs for interference terms. This is 
how the effect of quantum statistics enters the semiclassical worldline instanton approximation. The agreement with 
the WKB analysis of [28], which in turn agrees very well with exact results, is extremely good. 

We end with some brief comments about lessons from this work for the more general case, when the gauge field 
background represents a more complicated laser configuration, for example colliding short-pulse beams with spatial 
focussing. In this case, the old-style WKB approaches are not directly useful, but a possible semiclassical approach is 
provided by the worldline formalism. As emphasized already by Fock, Nambu, Feynman and Schwinger [3l [37l |42j |43] , 
the problem reduces to a problem of tunneling in four-dimensional Minkowski spacetime. As this work has shown, 
we need to consider complex classical paths, and furthermore the trajectories contain both instanton segments and 
interference segments. It would be interesting to investigate the applicability of various numerical and analytical 
methods that have been developed for multi-dimensional tunneling, in the context of non-relativistic two- and three- 
dimensional problems in chemical and molecular physics [3T1 US] . These methods also typically involve complex 
trajectories, in the spatial coordinates, and analytic continuation of the "time" coordinate to imaginary values. Such 
an extension would have to take into account the derivative coupling inherent to a gauge theory, whereas most previous 

-2 

results are for Hamiltonians of the form H = + V{x). (In the non-relativistic quantum mechanical context the 
analogue would correspond to tunneling in the presence of a magnetic field, as has been studied by Dykman |26|.) 
The QED extension would also have to take into account the relativistic causality features of Minkowski space, which 
result in the very different physical role played by electric and magnetic components of the background field. An 
interesting class of fields to investigate would be finite-plane-wave fields [T5], which involve both space and time, and 
which exhibit analogous quantum interference effects. 

We acknowledge support from the US DOE grant DE-FG02-92ER40716. 

VI. APPENDIX: AN IMPORTANT NUMERICAL TECHNICALITY 

For fields with localized structures, there arises an interesting numerical technicality due to the fact that in addition 
to turning points in the complex x° plane, there can also be poles of the function ^3(0;'^) in the complex x*^ plane. For 
example, for the field ^3(0;°) = E /uj tanh{ujx'^) these poles lie along the imaginary axis, at x'^ = {n + l/2)7r j/w, while 
for the field A3(x°) = E/u}/{l + (ujx^)^) there are precisely two poles, at x° = ±i/ijj. The locations of these poles do 
not depend on p^. However, as ps is varied, we have observed that beyond certain values of the turning points can 
move into positions where the classical trajectories, obtained by the numerical method outlined in Section IIIC, come 
close to the poles. In these cases, when the trajectory passes very close to the pole, the numerical integration of the 
equation of motion can lead to the trajectory "jumping" to another branch, so that it does not in fact continue to 
the expected final turning point. This applies to both the instanton and interference trajectories, depending on the 
location of the turning points relative to the poles. 

Similar numerical instabilities [although for slightly different reasons] have been observed in semiclassical studies 
of multi-dimensional quantum mechanical problems, where complex instantons are also required |34j . One possible 
solution is the prescription of Rubakov and collaborators: in fact, one only needs to minimize the imaginary part of 
the action in the path integral, which in turn means that the appropriate initial condition for instanton segments is 
to take just the imaginary part of the initial velocity to vanish. This amounts to a stationary phase approximation. 
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Then the real part of the initial velocity is a free parameter, and one should vary with respect to this parameter 
in order to find those trajectories that give a minimum imaginary part of the classical action. Implementing this 
prescription in our gauge theory case we find that in the situation where the classical trajectories do not approach 
the poles, our prescription produces completely equivalent results for the classical action on the trajectory segments. 
When the longitudinal momentum reaches a threshold value beyond which the trajectories "jump" to another branch, 
we can cure the situation by only fixing the imaginary part of the initial velocity to vanish. By adjusting the real 
part of the initial velocity, we can force the trajectory to go to the expected final turning point, and the correct 
minimum imaginary action is obtained by tuning the real part of the initial velocity to the threshold value where 
the jumping is first avoided. While this procedure can indeed be implemented, in these QED problems we found it 
also to be extremely delicate numerically. This is because one needs to compute the proper-time interval T of the 
trajectory segment very precisely in order to evaluate the classical action accurately. A small error in the period T 
can produce a significant error in the classical action, iWd, because the classical Lagrangian changes sign rapidly in 
the neighborhood of T. 

We have found another, more numerically stable, way to overcome the problem of classical trajectories approaching 
poles. The key observation is that we can relax the proper-time normalization condition (11), by taking advantage 
of the well-known reparametrization invariance of the path integral [311130]. Thus, for scalar QED, instead of ([6]), we 
have the more general expression: 



/■°° dT 

r— / 2?n$(n) — e- 



/„ n(u)du / ^4^ y g-iS 

Jx(0)=x(T) 



(46) 



where n represents the auxiliary einbein field, which is to be fixed with the aid of the <i>(n) "gauge-fixing" functional 
[5^ 140) . The action including the einbein field is: 



/o \ 2n clu au c 



(47) 



Previously we chose the " gauge- fixing" condition n — 1, but we are free to rescale proper-time by any factor. In 
particular, note that our choice of imaginary u for the instanton segments can be thought of as taking n = i. 




Im(xo) 
10 ■ 
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Re(xo) 
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FIG. 9: Producing closed trajectories by tuning the einbein n, as explained in the text. These plots are for the two-turning- 



point-pair gauge field in (41 1. The first plot shows the situation for Eq — 0.1, lj = 0.1, and pa — 1. The interference trajectories 
do not connect the different turning point pairs, due to the appearance of poles at = ±i/uj = ±10i. The second plot shows 
the result of choosing n — 1 ± 0.74i and integrating the modified equations ( |49| |. Now the poles are avoided and a closed 
classical trajectory results. The action |50| evaluated on this trajectory gives the correct semiclassical approximation to the 
particle spectrum, as plotted in Figure!^ 

Including this einbein factor, the classical equations of motion, together with the constant of integration, read: 

^=n-F^^x)^ , ^^^n'm'c^ (48) 
dw^ c du du du 



Therefore, instead of numerically integrating the equations of motion ( [3l] [32] [33| and ([34] [35j [36] ), we can integrate 
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the following coupled equations, with associated initial conditions: 



dv? 
dv? 

("initial) 
("initial) 

d^ 

du 

dX3 

du 



-n- — doAsix'') 
c du 



-n — dMx'') 
c du 

„o 



a turning point solution of ( 30 ) 



(49) 



The choice n = 1 gives the interference segments from before, while the choice n = i gives the instanton segments 
as before. Now, to find the instanton segments we choose n = i{l + ib), and for the interference segments we choose 
n = (1 + ib), for some real b. We can tune b so that the trajectories avoid the poles. Finally, given such a trajectory, 
we compute the associated classical action as 



W[x°iuy,-m^c^] = 



dx° 
du 



du 



(50) 



As an illustration of this procedure, consider the field (41) with two sets of complex conjugate turning points. In the 
first plot in Figure [9] we see that for P3 — 3 the interference trajectories do not connect the different pairs of turning 
points, as expected, but go off to infinity. This is with the choice n = 1. We can cure this by tuning n to take the value 
n = 1 ± 0.74i, with the resulting trajectories shown in the second plot of Figure |9] This choice of n is chosen so that 
the classical trajectories pass safely by the poles that occur at x° = ±i/w [in the plot, l/ui = 10]. They are patched 
together smoothly on the imaginary x'^ axis. Now we have the desired instanton segments, and after evaluating the 



classical action ( |50[ ), we obtain the correct semiclassical contribution to the particle number momentum spectrum, 
as shown in Figure [8j Note that for this — 3, and n — i, the instanton segments do indeed connect the complex 
conjugate turning points. We have found that by suitable tuning of the einbein factor n we are always able to avoid 
the poles and produce closed classical trajectories that avoid the poles. 
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